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We report the detection of highly-ionized gas at z ~ seen in resonant UV 
and X-ray absorption lines toward the z = 0.03 blazar Mkn 421. A total of 13 
X-ray and 3 UV lines were measured (or upper limits derived), including three 



lines in the O VII K-series and Ka transitions from neon, carbon, and nitro- 

00 ' 

gen. From the three VII lines we derive a 2a Doppler parameter constraint 
of 24 < b < 55kms _1 . The FUSE spectrum shows strong Galactic low-velocity 
O VI 1032 A absorption as well as a possible weak O VI high-velocity component 
(HVC). The Doppler parameter of the low-velocity O VI measured with FUSE 
is ~ 3(7 higher than that derived from the O VII line ratios, indicating that 
CLf the O VII and Galactic O VI arise in different phases. This velocity dispersion, 

along with limits on the gas temperature and density from the X-ray line ratios 
(assuming a single phase with collisional ionization equilibrium plus photoion- 
ization) are all consistent with an extragalactic absorber. However, the O VII 
Doppler parameter is inconsistent with the high temperature required to produce 
the observed O VIhvc/O VII ratio, implying that the HVC is probably not re- 
lated to the O VII. In addition, the O VI k a line detected by Chandra implies a 



department of Astronomy, The Ohio State University, 140 West 18th Avenue, Columbus, OH 43210, 
USA 

2 Harvard- Smithsonian Center for Astrophysics, Cambridge, MA, USA 

3 Instituto de Astronomia Universidad Autonomica de Mexico, Apartado Postal 70-264, Ciudad Univer- 
sitaria, Mexico, D.F., CP 04510, Mexico 

4 University of California-Berkeley, Berkeley, CA, USA 

5 Chandra Fellow 

6 Osservatorio Astronomico di Roma, Monteporzio, RM, Italy 

7 Department of Astronomy, University of Massachusetts, Amherst, MA 01003, USA 



- 2 - 



column density ~ 4 times higher than the 1032 A absorption. Although an extra- 
galactic absorber is fully consistent with the measured column density ratios, a 
Galactic origin cannot be ruled out given the uncertainties in the available data. 

Subject headings: Galaxy: halo — intergalactic medium — quasars: absorption 
lines 



1. Introduction 

At both low and high redshifts, the intergalactic medium (IGM) is thought to be the 
dominant reservoir of baryonic matter in the universe. Essentially all of the baryons at z ~ 2 
are accounted for in the "forest" of H I Lyman-a absorption lines seen against background 
quasars (Weinberg et al. 1997). In the local universe, numerical simulations predict that most 
of the IGM is in the form of a shock-heated "warm-hot IGM" (WHIM) with T ~ 10 5 - 10 7 K 
and densities n H ~ 10~ 6 - l(T 5 cur 3 (Cen & Ostriker 1999; Dave et al. 2001). Hydrogen is 
mostly ionized at these temperatures, rendering it difficult or impossible to detect. However, 
H-, He-, and Li-like ions of heavier elements can exist under these conditions, producing a 
"forest" of \ow-z intervening UV/X-ray absorption lines in quasar spectra (Hellsten et al. 
1998; Perna & Loeb 1998). In particular, O VI, O VII, and O VIII are expected to produce 
the strongest lines (depending on the temperature and ionizing background). 

The required sensitivity to detect such absorption systems is now available with the 
advent of such facilities as the Chandra X-ray observatory, XMM-Newton, Hubble Space 
Telescope, and the Far Ultraviolet Spectroscopic Explorer (FUSE). High signal-to-noise ob- 
servations along four lines of sight reveal such intervening absorption systems with tempera- 
tures and densities consistent with the WHIM (Nicastro et al. 2005a,b,c; Mathur et al. 2003; 
Fang et al. 2002, Fang et al., in prep). Additionally, the total baryonic mass density in these 
detected filaments is consistent with the (2.5±0.4)% deficiency between CMB measurements 
(Bennett et al. 2003) and previous density measurements of stars and gas (see also Nicastro 
et al. 2005a, and references therein). 

Since most galaxies are expected to trace the same cosmic overdensities as the web 
of IGM filaments, it would be no surprise if the Galaxy itself resided in such a filament. 
Indeed, X-ray spectra of several quasars show likely z = O VII absorption, e.g. PKS 2155- 
304 (Nicastro et al. 2002), 3C 273 (Fang et al. 2003), and Mkn 279 (Williams, Rodriguez et 
al., in preparation). The upper limit found for O VII emission implies that O VII absorption 
systems most likely arise in an extragalactic medium (Rasmussen et al. 2003); however, the 
analysis is complicated by z = O VI absorption observed with FUSE along many quasar 
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lines of sight (Wakker et al. 2003). The strongest O VI absorption component, typically 
within t>LSR = ±100 km s" 1 (hereafter denoted O VIlv), is due to gas in the halo and thick 
disk of the Galaxy (Savage et al. 2003). Of the 102 sightlines studied, 84 show at least one 
additional O VI high-velocity component (HVC), typically weaker than the O VIlv- 

The origin and location of the O VI HVCs is controversial: while some may be spatially 
and kinematically associated with known Galactic structure (such as the Magellanic stream 
and H I high-velocity clouds, though this association is in many cases uncertain), others 
appear to be completely isolated. Sembach (2003) summarizes the evidence in favor of a 
Galactic origin for the O VI HVCs. In such a scenario, the absorption could arise at the 
interface between infalling cooler gas clouds and very hot (T > 10 6 K) gas in an extended 
Galactic corona. The observed gas stripping and orbital deflection of the Magellanic clouds 
also lends credence to such a low-density corona. Furthermore, lower-ionization species 
such as Si IV and C IV are seen toward some sightlines. Thus, there is some evidence for a 
collisionally-ionized, high-density Galactic medium. 

On the other hand, there is also evidence that the isolated O VI HVCs trace the Local 
Group medium and thus are associated with the z = O VII and O VIII X-ray absorption. 
The most compelling evidence for this interpretation is the velocity segregation of the O VI 
HVCs: Nicastro et al. (2003) report that the mean velocity vector of the O VI HVCs is 
highest in the local standard of rest and lowest in the Local Group rest frame, implying 
an extragalactic, Local Group origin. The presence of lower-ionization lines does not rule 
out a WHIM origin for the O VI HVCs since the medium is likely inhomogeneous and may 
thus consist of phases of different temperatures and densities (such as the intervening WHIM 
absorbers found by Danforth & Shull 2005). Furthermore, isolated high-velocity O VI (which 
is not associated with O I and other low-ionization lines) has not been detected toward high 
latitude Galactic sources out to lOkpc (Danforth et al. 2002; Howk et al. 2003; Zsargo et 
al. 2003). The low-ionization line profiles often differ substantially from the O VI HVCs, 
implying that they may arise in different components (Collins et al. 2005). 

The key to understanding the O VI absorption is higher-quality X-ray data, since the 
O VII and O VIII along high-latitude quasar sightlines are the most likely to trace the WHIM 
(i.e., least subject to possible contamination by hot Galactic gas). X-ray measurements by 
themselves provide great insight into the properties of the local WHIM, including its mass 
and extent, thereby yielding valuable constraints for models of galaxy formation and the 
Local Group (e.g., Kravtsov et al. 2002). However, knowing whether the O VI HVCs are 
also associated with the WHIM provides even stronger diagnostics. Unfortunately, existing 
X-ray data have not been of sufficiently high quality to resolve this ambiguity: even at large 
column densities (iVj ~ 10 15 cm~ 2 ), X-ray lines are quite weak, and several thousand counts 
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per resolution element in a typical Chandra grating spectrum are required to accurately 
measure these lines. 

Through our program of observing nearby blazars in outburst phases, we have obtained 
high-quality Chandra and FUSE spectra of Mkn 421, sufficient to study in detail the local 
WHIM (and Galactic halo/thick disk) absorption. Here we report on these observations, 
and the inferred properties of the local absorption. 

2. Observations and Data Preparation 

2.1. Chandra 

A full description of the Chandra observations, data reduction, and continuum fitting can 
be found in Nicastro et al. (2005a); a brief summary follows. Mkn 421 was observed during 
two exceptionally high outburst phases for 100 ks each as part of our Chandra-hOA observing 
program: one at /o.5-2kcV — 1-2 x 10~ 9 ergs -1 cm~ 2 with the Low Energy Transmission 
Grating (LETG) combined with the Advanced CCD Imaging Spectrometer-Spectroscopic 
(ACIS-S; Garmire et al. 2003) array, and another at /o.5-2keV = 0.8 x 10~ 9 ergs -1 cm~ 2 
with the High Resolution Camera-Spectroscopic (HRC-S; Murray & Chappell 1985) array 
and LETG. Each of these observations contains ~ 2500 counts per resolution element at 
21.6 A. Additionally, another short observation of Mkn 421 was taken with HRC/LETG (29 
May 2004), providing another 170 counts per resolution element. These three spectra were 
combined over the 10-60 A range to improve the signal-to-noise ratio (S/N~ 55 at 21 A with 
0.0125 A binning). The final coadded spectrum of Mkn 421 is one of the best ever taken 
with Chandra: it contains over 10 6 total counts with ~ 6000 counts per resolution element 
at 21.6 A, providing a 3a detection threshold of W\ ~ 2 mA (AT OVII = 8 x 10 14 cm~ 2 for an 
unsaturated line). 

Effective area files (ARFs) for each observation were built using CIAO 1 v3.0.2 and 
CALDB 2 v2.2.6. Those pertaining to the ACIS/LETG observations were corrected for the 
ACIS quantum efficiency degradation 3 (Marshall et al. 2003). For the HRC/LETG obser- 
vations, the standard ARFs were used. Each ARF was then convolved with the relevant 
standard redistribution matrix file (RMF), and the convolved RMFs were weighted by ex- 

1 http://cxc. harvard.edu/ciao/ 
2 http://cxc. harvard.edu/caldb/ 

3 See also http://cxc.harvard.edu/ciao/why/acisqedeg.html 
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posure time, rebinned to the same energy scale, and averaged to provide a response file for 
the coadded spectrum. 

Using the CIAO fitting package Sherpa 4 , we initially modeled the continuum of Mkn 421 
as a simple power law and a Galactic absorbing column density of N H = 1.4 x 10 20 cm~ 2 
(Dickey & Lockman 1990), excluding the 48-57 A HRC chip gap region. Metal abundances 
for the Galactic gas were then artificially adjusted to provide a better fit around the O I 
and C I K-edges near 23 A and 43 A respectively. This is not intended to represent actual 
changes to the absorber composition, but rather to correct uncertainties in the instrument 
calibration. These adjustments affect the continuum mostly near the carbon, oxygen, and 
neon edges, but individual narrow absorption lines are unaffected. After this fit there were 
still some systematic uncertainties in the best-fit continuum model; these were corrected 
with broad (FWHM = 0.15 - 5 A) Gaussian emission and absorption components until the 
modeled continuum appeared to match the data upon inspection. Indeed, the residuals of the 
spectrum to the final continuum model have a nearly Gaussian distribution, with a negative 
tail indicating the presence of narrow absorption lines (see Nicastro et al. 2005a, Figure 8). 



2.2. FUSE 

Mkn 421 was also observed with FUSE as part of our Director's Discretionary Time 
observing program on 19-21 January 2003 for a total of 62.8 ks. An additional 21.8 ks 
observation from 1 December 2000 was also available in the archive. We used the time- 
tagged, calibrated data from only the LiFlA detector, since inclusion of the LiF2B data 
provides a small (~ 20%) increase in S/N but degrades the overall spectral resolution. 5 
These two observing programs consist of four observations, which in turn contain a total of 29 
individual exposures. The wavelength scales of each observation's constituent exposures were 
cross-correlated and shifted (typically by 1-2 pixels) to account for slight uncertainties in the 
wavelength calibration. The exposures for each observation were checked for consistency and 
coadded, weighted by exposure time. The resulting four spectra were then cross-correlated 
against each other, coadded (with a ~ 10% downward shift in flux applied to the 2000 
observation due to source variability), and binned by 5 pixels (0.034 A, or one-half of the 
nominal 20kms _1 resolution) providing a S/N of 17 near 1032 A. 

To check the absolute wavelength calibration we followed the method of Wakker et al. 



4 http://cxc. harvard.edu/sherpa/ 

5 See the FUSE Data Analysis Cookbook vl.0, http://fuse.pha.jhu.edu/analysis/analysis.html 
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(2003), using their 4-component fit to the Murphy et al. (1996) Green Bank H 1-21 cm 
data as a velocity reference. They find four main components of H I emission with an Nr- 
weighted average velocity of -31.7 km s -1 . In the FUSE spectrum, the Si II A1020.699 A 
and Ar I A1048.220A lines are expected to trace the same gas as the H I emission. Each 
UV line was fit with two Gaussian components in Sherpa, giving average velocity offsets of 
—30.9 and — 34. 9 km s^ 1 respectively. These agree well with the H I data, though the slight 
difference between the Ar I and Si II measurements suggest at least a ~ 4kms _1 intrinsic 
wavelength uncertainty. 



3. Line Measurements 

To find and identify narrow absorption lines in the Chandra spectrum of Mkn 421, we 
visually inspected small (2-5 A) regions of the spectrum, beginning around the rest wave- 
length of O vn Ka (21.602 A) since this tends to be the strongest z = X-ray absorption 
line (e.g. Nicastro et al. 2002; Fang et al. 2003; Chen et al. 2003). Three O VII^ a (Figure 1) 
absorption features were found: one at z — 0, one at z — 0.011, and one at z — 0.027 (with 
typical redshift errors of 0.001). There is also a strong feature which is ~ 3a from the O Y Ka 
rest wavelength, but is more likely O VII^ a at v ~ +900 km s -1 relative to the blazar. A 
close pair of lines consistent with Lja at this velocity has been observed (Shull et al. 1996; 
Penton et al. 2000), so this may be indicative of an inflow to Mkn 421 or uncertainty in 
the blazar redshift (based on rather old spectrophotometric measurements by Margon et 
al. 1978). A weak O VI^ a line is seen at the rest wavelength of 22.02 A. Other regions of 
the spectrum were then searched for lines corresponding to these systems, with particular 
emphasis paid to strong transitions of the most abundant elements (C, N, O, and Ne). All 
in all there were 13 lines marginally or strongly detected at z ~ (including the N VII, O V, 
and Ar XV upper limits), 3 at z = 0.011, and 7 at z = 0.027. The latter two systems are the 
subject of other papers (Nicastro et al. 2004a,b) and thus will not be discussed further here. 

These 13 z = X-ray lines were fitted in Sherpa with narrow Gaussian features su- 
perposed on the fitted continuum described in §2.1 (see Figure 1). We are excluding the 
strong O I (23.51 A) line since it arises in the neutral ISM and is not of interest here, as 
well as the O2 (23.34 A) absorption since it coincides with a strong instrumental feature 
and cannot be accurately measured. Due to the FWHM = 0.04 A (~ 600 km s" 1 ) LETG 
resolution the lines are all unresolved, so only the position and equivalent width of each line 
are measured. Errors are calculated using the "projection" command in Sherpa, allowing the 
overall continuum normalization to vary along with all parameters for each line. The result- 
ing line parameter estimates are presented in Table 1. The ~ 0.02 A systematic wavelength 
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uncertainty of the LETG 6 is in most cases larger than the statistical uncertainty of the line 
centroid; thus, Table 1 lists whichever is greater. Additionally, a meaningful upper error bar 
on the C VI equivalent width could not be calculated with Sherpa. In this case, the FWHM 
was frozen at the instrumental resolution and the error was recalculated; a visual inspection 
confirms the new limit to be more reasonable. Upper limits for the O V, N VII, Ar XV, and 
Ne X lines were calculated with both the position and FWHM frozen. 

The FUSE spectrum (Figure 2) shows a strong, broad low-velocity O VI 1032 A ab- 
sorption line at z ~ due to gas in the Galactic thick disk and halo (Savage et al. 2003). 
An asymmetric wing on the red side of this line is evident, possibly a kinematically distinct 
HVC. We fitted the O VI 1032 A line in Sherpa using a constant local continuum (in a ±2 A 
window) and two Gaussian absorption components: one for the t)~0O VIlv line, and one 
at v ~ 100 km s -1 for the HVC. No H 2 contamination is seen at the O VI 1032 A wavelength 
when absorption templates are fit to the other H 2 lines seen in the spectrum. The 1037 A line 
is somewhat blended with a single H 2 absorption line; this is taken into account with another 
narrow Gaussian. From this fit, we find equivalent widths of 18.6 ± 5.6 mA for the 1032 A 
HVC and 270.7 ± 7.9 mA for the Galactic component. The best-fit model for the HVC is 
fairly robust and not sensitive to variations in the initial parameters; however, the derived 
equivalent width of 18.5 ± 5.6 mA is lower than the 37 ± 11 ± 29 mA (errors are statistical 
and systematic, respectively) measured by Wakker et al. (2003) in the initial 21.8 ks obser- 
vation. They employed a direct integration method which may not have taken into account 
the substantial blending of the Galactic O VI L v with the HVC. Our total O VI equivalent 
width (LV+HVC = 279 ± 10 mA) is in good agreement with their value of 285 ± 20 mA. 

Deb lending the O VI 1037 A line is less certain due to the presence of adjacent Galactic 
C II and H 2 absorption. A flat continuum was again employed from 1035 — 1040 A and 
single Gaussian components were used to fit the C II, O VI, and H 2 absorption lines. The 
HVC on the 1032 A O VI line should also appear at ~ 1038 A with W x = 0.50 x W A (1032 A). 
Although this component is too weak to be detected directly, it could cause the measurement 
of the Galactic LV-0 VI 1037 A line to be systematically high. Another absorption Gaussian 
with one-half of the 1032 A HVC equivalent width (and with the same FWHM and velocity 
offset) was included in the 1037 A O VI line fit to account for this. Table 1 lists the measured 
properties of the O VIlv and HVC absorption lines. 



6 http://cxc. harvard.edu/cal/ 
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4. Absorption Line Diagnostics 

4.1. Doppler Parameters 

To convert the measured equivalent widths to ionic column densities, we calculated 
curves of growth for each absorption line over a grid of Doppler parameters {b = 10 — 
100 km s -1 ) and column densities (logiV^/cm -2 = 12.0 — 18.0), assuming a Voigt line profile. 
Since the X-ray lines are unresolved, b cannot be measured directly. It can, however, be 
inferred from the relative strengths of the three measured O VII K-series lines. These lines 
are produced by the same ionic species, so in an unsaturated medium W\ oc fi u \ 2 where fi u 
is the oscillator strength. The expected equivalent width ratio of O VII Kf3 to Ka is then 
W\(KP)/W\(Ka) = 0.156, so the measured value of 0.49 ±0.09 indicates that the Ka line is 
saturated. On the other hand, the measured O VII K7/K/5 ratio is 0.43 ±0.16, in agreement 
with the predicted (unsaturated) value of 0.34. 

These line ratios by themselves are insufficient to determine the physical state of the 
O VH-absorbing medium since b and iVovii are degenerate: the Ka line saturation could be 
due to high column density, low b, or a combination of both. However, given an absorption 
line with a measured equivalent width and known fi u \ 2 value, the inferred column density 
as a function of the Doppler parameter can be calculated. The measured equivalent width 
(and errors) for each transition thus defines a region in the iVovii — b plane. Since the actual 
value of iVovii is fixed, b and Aovii can be determined by the region over which the contours 
"overlap;" i.e. the range of Doppler parameters for which the different transitions provide 
consistent iVovii measurements. 

Figure 3 shows such la contours for the three measured O VII transitions. As expected, 
the inferred iVovii is nearly constant in the unsaturated regime (large b), and rises sharply 
at low b as the lines begin to saturate. At each value of b, the differences A(logA r Q/ g) = 
log(NK a ) — loglNxp) and A(logiV Q7 ) = log(iV^ Q ) — log(iV^ 7 ) were calculated, along with 
the errors on each A(logiV). The quantity A(logA r Qj a) is consistent with zero at the la and 
2a levels for 15 < b < 46kms _1 and 13 < b < 55kms _1 respectively, while A(logiV Q , 7 ) 
provides limits of 31 < b < 50kms _1 and 24 < b < 76kms _1 respectively. Since A(logiV Q , 7 ) 
provides a more stringent lower limit on b while A(log N a/3 ) better constrains the upper limit, 
we thus assume a la range of 31 < b < 46kms~ 1 , and a 2a range of 24 < b < 55kms _1 . 
It should be noted that Figure 3 also shows some overlap between the K a and K 1 at b < 
12kms _1 ; however, this solution is unlikely given the lower limit provided by the Kp line. 
Moreover, b = ^kms" 1 implies a maximum temperature (assuming purely thermal motion) 
of T max = 1.3 x 10 5 K; such a low temperature is unlikely to produce the observed strong 
high-ionization lines. 
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A similar analysis is not as effective when applied to the strong VIlv UV doublet 
(from the thick disk), since these lines are only slightly saturated: the measured W\ ratio 
is 0.61 ± 0.04, compared to the expected unsaturated value of 0.50. When the inferred 
iVovi is calculated as a function of b for both lines of the O VIlv doublet, the predicted 
iVovi values are consistent over b = 34 - 112 km s" 1 (at the 2a level; see Figure 4). Since 
the O VIlv 1032 A line is fully resolved by FUSE (~ 15 resolution elements across the 
line profile) and relatively unblended, its Doppler parameter can be estimated much more 
accurately using the measured line width and strength. In an unsaturated absorption line, 
FWHM = 2(ln2) 1/2 6; however, the measured FWHM increases if the line is saturated. 
We compensated for this by calculating Voigt profile FWHMs on a grid of iVovi and b, 
and determining the region consistent with the O VIlv 1032 A FWHM measurement of 
152±7kms~ 1 (or b = 91 ±4kms~ 1 assuming no saturation). 

When the FWHM-derived contour is overlaid on the Aqvi — b contour inferred from the 
equivalent width measurement of the LV-0 VI 1032 A line, the two regions overlap nearly 
orthogonally (Figure 4) leading to a constraint of 6(0 VIlv)= 80.6 ±4.2 km s" 1 . This is ~ 2a 
lower than the unsaturated FWHM, once again confirming that the O VIlv is only weakly 
saturated. At this b the inferred column densities from the two lines of the O VIlv doublet 
differ by 1.3c but this is only a minor discrepancy and likely due to errors introduced by the 
blending of the 1037 A line; thus, we will only consider results from the more reliable 1032 A 
line measurement. However, at no value of the Doppler parameter do the 1032 A, 1037 A, 
and O Vlxa lines all produce a consistent iVovi measurement; in fact, the O Vlxa column 
density is a factor of ~ 4 higher than that inferred from the UV data. This discrepancy is 
discussed further in §5.1.1. 



4.2. Column Densities 

The Doppler parameters measured for the O VII (31kms _1 < b < 46kms _1 ) and 
O VIlv (b = 80.6 ± 4.2 km s -1 ) absorption are inconsistent at the ~ 3a level, indicating the 
presence of at least two distinct components: the Galactic thick-disk gas traced by broad 
t)~0O VIlv absorption, and another lower-6 phase, possibly of extragalactic origin, traced 
by the O VII absorption lines. It cannot be assumed a priori that any given line (other than 
those used to determine b) originates in one phase or another; moreover, the uncertainty in 
the calculated column density depends not only on the equivalent width error but also the 
error in b. To take this into account, for each ion the derived column density logA^ and 
its ±lcr limits were averaged over the ±ler ranges of both measured Doppler parameters. 
As it turns out, the choice of b does not make a significant difference since all other lines 
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(besides the O VII and O VIlv absorption) are essentially unsaturated; i.e., the difference in 
Ni calculated with the O VIlv and O VII Doppler parameters is small compared to the la 
error on the equivalent width measurements. Even so, to avoid possible systematic errors, 
we assumed b = 80.6 ±4.2 km s -1 for those lines likely to originate in the Galactic thick disk 
(O VI LV , O V, and C V), and b = 31 - 46kms~ 1 for all other species. The derived ionic 
column densities are listed in Table 1 



4.3. Temperature and Density Constraints 

At densities such as those found in the Galactic interstellar medium (ISM; n e ~ 1 cm -3 ), 
photoionization is unimportant because thermal collisions are by far the dominant ioniza- 
tion source. This is also the case for very high temperatures (T > 10 7 K) even at low 
densities, since the collisional rate is greater than the photoionization rate. However, at the 
low densities typically found in the intergalactic medium (n e = 10~ 6 — 10~ 4 cm -3 ), photoion- 
ization from the diffuse UV/X-ray background begins to play a greater role by enriching the 
abundances of highly-ionized elements at typical WHIM temperatures (logT(K) ~ 5 — 7) 
relative to those expected from pure collisional ionization (Nicastro et al. 2002; Mathur et 
al. 2003). It is thus imperative that the ionizing background be taken into account in order 
to accurately predict ionic abundances in the WHIM. 

Version 90.04 of the ionization balance code Cloudy (Ferland 1996) was used to compute 
collisional- plus photoionization hybrid models for the absorbing medium. Relative ionic 
abundances were computed over a grid of logT(K) = 4.5 — 7.4 and logn e (cm~ 3 ) = —7 — 0, 
with a step size of 0.1 dex in both logn e and logT. Initially, a rigid scaling of [Z/H] = — 1 
for all metals was assumed. For the ionizing background we employed the Sternberg et al. 
(2002) fit to the metagalactic radiation field: 

{ Juo (-Y 13 K - < 4 

Ju=\ W f ,-0.46 •* (1) 

I 2-512 x 10- 2 J,o (£) £ > 4 

where here J u0 = 2 x 10~ 23 ergss _1 cm -2 Hz -1 sr -1 and u = 13.6 eV. The total flux of 
ionizing photons is then given by f 1 = An J™{J v /hv)dv = 1.3 x 10 4 photons s _1 cm -2 , and 
the ionization parameter is \ogU = log(/ 7 /c) — logn e = —6.36 — log(n e ) where n e is the 
electron density in cm -3 . 

Using the ionic abundances calculated with Cloudy, we derived expected abundance 
ratios for all observed ions at each point in the logn e — logT plane. Since any given density 
and temperature uniquely determines a set of abundance ratios (N a /Nb for all ions a and 
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fe), the problem can be inverted: any value of N a /N b defines a curve in the logn e — logT 
plane, i.e. a set of temperatures and densities which can produce the measured ratio. When 
the errors on jV a /iVj, are taken into account, the curves become contours, and the overlap 
between two or more contours defines the temperatures and densities for which the measured 
ratios are consistent. This is analogous to the method used in §4.1 to determine Doppler 
parameters for O VIlv and O VII. 

The most powerful diagnostics are those using ratios between different ions of the same 
element, since these ratios are independent of the relative metal abundances. Unfortu- 
nately the N VII/N VI and Ne X/Ne IX upper limits are not stringent enough to place 
meaningful constraints on the temperature and density. Figure 5 shows the logn e — logT 
contours for ratios between the X-ray O VI^ Q , O VII, and O VIII lines as well as the 
O VIhvc/O VII ratio. The X-ray line ratios are inconsistent with a high-density (n e > 
10~ 3 cm~ 3 ), high-temperature (logT > 6.2) medium, and instead converge on a partially pho- 
toionized plasma with n e = 10~ 4 ' 7 — 10~ 3 ' 9 cm~ 3 (from the overlap between the O VIxa/O VII 
and O VIII/O VII contours) and T = 10 5 - 5 ~ 5 - 7 K (from the limits provided by O Vl Ka /0 VII 
in this density range). These ranges of temperatures and densities are in line with those 
expected from WHIM gas (Dave et al. 2001). Of course, this is all contingent on the O Ylxa 
line being a reliable tracer of iVovi; this caveat is discussed in detail in §5. 

On the other hand, the O VIhvc/O VII ratio overpredicts the temperature by at least an 
order of magnitude for all values of logn e — in order to be consistent with the O VIII/O VII 
ratio, the O VIhvc/O VII ratio would need to be stronger by a factor of ~ 2.5 (or ~ 3. 5c). It 
is possible that the HVC is not a physically distinct component, but is instead the result of 
some systematic error (such as fixed pattern noise or an unexpected anomaly in the Galactic 
O VI L v velocity distribution). In this case, the O VI associated with the O VII and O VIII 
may be completely blended with the thick-disk O VIlv and thus unmeasurable. Consistency 
with the O VIII/O VII ratio (in the collisional ionization regime) requires log(iV viA / Vovii) ~ 
—2.5, or roughly 20% of the Galactic UV O VI absorption. Although it appears that the 
O VI HVC as measured cannot originate in the same medium as the O VII absorption, we 
suspect that additional atomic physics may be at work here and could in principle reconcile 
this disagreement (see §5.1.1). 

While the O VIr-q/O VII and O VIII/O VII ratios provide strong constraints, it is 
important to consider other ion ratios as well (particularly since the O VI^ a and 1032 A 
O VI column densities disagree). Figure 6 shows the logn e — logT contours for several 
different ion ratios, all calculated relative to O VII since the error on iVovn is small. With a 
rigid metallicity shift relative to Solar, the Ne IX/O VII O VIII/O VII, and O VI^ a /0 VII 
ratios are not all consistent with each other for any combination of temperature and density; 
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however, the consistency can be improved with adjustments to the [Ne/O] ratio (see §5.1.3). 
Both the C VI/O VII and N VII/O VII measurements are consistent with a high- or low- 
density medium at solar abundances. 

Limits on the temperature of the Galactic thick-disk absorption can be derived in a 
similar fashion, although it is not the primary focus of this work and there are far fewer 
measured lines to work with. The most accurately-measured line is the O VIlv; additionally, 
C V and N VI X-ray lines are measured, and upper limits have been determined for O V 
and N VII. Figure 7 shows the temperature constraints derived for this Galactic absorption, 
assuming pure collisional ionization. The O V/O VIlv and N VII/N VI upper limits provide 
metallicity-independent constraints of logT > 5.39 and logT < 6.64 respectively. A more 
stringent upper limit on temperature of logT < 6.03 is provided by the N VII/O VI L v 
ratio, but this is somewhat dependent on [N/O]. Within this range the C V/O VI L v ratio 
provides an even stricter limit of 5.3 < logT < 5.7, but again this depends on [C/O]. At 
these temperatures the expected O VII column density is at most an order of magnitude less 
than measured; thus, the O VIlv, O VII, and O VIII cannot all originate in the same phase 
assuming pure collisional ionization (see also Mathur et al. 2003). 

5. Discussion 

Our Chandra and FUSE observations have provided a wealth of data on absorption 
near the Galaxy, constraining the temperature and density tightly (logT(K) = 5.5 — 5.7 and 
n e = 10~ 4 ' 7 — 10~ 3 ' 9 cm~ 3 when the O Vl^a measurement is included), which are conditions 
suggestive of the local group intergalactic medium and require supersolar [Ne/O]. Here we 
first examine the assumptions that have led us to these results (§5.1), and then we discuss 
their implications for the location of the absorbing gas (§5.2), subject to these caveats. 

5.1. Potential Caveats 

5.1.1. The O VI Discrepancy 

The interpretation of the UV and X-ray data are particularly important, since (as 
Figure 5 shows) the combined O VIII/O VII and O VI/O VII ratios can provide tight 
constraints on the absorber temperature and density simultaneously (see also Figure 5 in 
Mathur et al. 2003). However, in this case the O VI column density inferred from the O YIro 
is a factor of ~ 4 higher than the combined 1032 A low- and high-velocity components. Since 
both the X-ray and UV transitions trace the same atomic state, the inferred column densities 
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should match. A similar disagreement has been seen in intrinsic AGN absorption systems 
(see Krongold et al. 2003; Arav et al. 2003); however, in these cases it is typically attributed 
to saturation or a velocity-dependent covering factor, neither of which is relevant to this 
z ~ absorption. 

On the other hand, our O Wxa measurement provides a test for these attributions; the 
local absorption, after all, is likely a dramatically different physical system than an AGN 
outflow, yet the same conflict arises. A macroscopic explanation does not adequately describe 
how this discrepancy is seen in both physical systems, so the actual reason may lie in the 
atomic physics of highly ionized plasmas. For example, some fraction of the O VI may be 
excited through collisions or recombination from O VII, and thus unable to produce 1032 A 
absorption while still absorbing O Vl^a photons. While a scenario that produces significant 
depopulation of the O VI ground state is difficult to envisage in such a low density plasma, 
we are investigating further the statistical equilibrium of O VI including photoexcitation 
and recombination in order to study such effects in more detail. However, it should be 
emphasized that this is not an isolated case so there must be a physical explanation for the 
O VI discrepancy, and the resolution of this paradox is crucial to our understanding of O VI 
UV and X-ray absorption and how it relates to the O VII. 

There is also the possibility that the line was misidentified as O Vlxa, and is actually 
another intervening O VII absorption line at z = 0.0195. This latter explanation is unlikely 
since no other absorption lines at this redshift are seen in the FUSE or Chandra spectra; 
additionally, this would require the line to fall exactly on the O VI rest wavelength, which 
seems like an improbable coincidence. Another possibility is that the theoretical oscillator 
strength of the O Vlxa transition is incorrect, but the value given in Pradhan (2000) would 
need to be low by a factor of ~ 2 — 4, in sharp contrast to the successful calculations of f tu for 
inner shell transitions in other ions in the same paper. Nevertheless, due to the discrepancy 
between the UV and X-ray O VI column density measurements, we present both possibilities: 
either (a) the O Vlxa line measures iVovi, or (b) it does not and is thus ignored. 

5.1.2. Absorption Components 

The Doppler parameter measurements indicate the existence of two distinct components 
along the Mkn 421 line of sight: one seen in the thick-disk O VIlv 1032 A absorption with 
^lv — 80.6 ± 4.2kms _1 , and the O VII absorber with foovn — 31 — 46kms~ 1 (la limits). 
The O VI HVC may represent a third phase (if case (b) above is correct) with &hvc = 
35^Qkms _1 (from the FWHM measurement). This agrees surprisingly well with the O VII 
b measurement, and is consistent with numerical simulations of the nearby IGM (Kravtsov 
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et al. 2002). However, the extremely low O VIhvc/O VII ratio requires a temperature much 
higher than the upper limit provided by feovn- 111 order for the HVC to trace the same gas as 
O VII (case a), then, the aforementioned atomic physics effects would need to be suppressing 
O VI HVC absorption and not the O VIlv line. Sembach et al. (2003) list mean Doppler 
parameters for a variety of HVCs, both Galactic and probable Local Group; unfortunately, 
the dispersion in these values and the errors on fcovn an d &hvc measured here are both too 
large to associate the components presented here to one of their classifications. 

It is also important to note that our analysis assumes a single phase origin for the 
included X-ray lines. This assumption is consistent with the data, given the good agreement 
between the three O VII lines in the calculated ranges of b and iVovi (Figure 3). Even 
so, if any of the ionic species arises in more than one phase along the line of sight, our 
results could be affected. For example, a Galactic, purely collisionally ionized medium can 
in principle reproduce the observed relative abundances of O VIII, O VII, O Vlxa, and 
Ne IX if several unresolved components are invoked to explain this absorption. However, 
the simplest explanation (a single-phase, low-density, partially photoionized extragalactic 
absorber) is fully consistent with all of these line measurements, and the similarity of the 
derived absorber properties to expectations for the local WHIM lend further support to this 
model (see §5.2). 



5.1.3. Abundances 

Due to the uncertainty associated with the O Vlxa absorption, metal abundances rel- 
ative to oxygen play a particularly important role in this analysis. By adjusting the metal 
abundances of the absorbing gas, the consistency of the solutions can in principle be im- 
proved with and without the O Vlxa measurement. Although the logn e — logT contour 
plots are useful for finding solutions, they cannot easily be used to determine the effects 
of changes in elemental abundances; thus, Figure 8 shows Ni/N y U as a function of logT 
for both the logn e = —3.9 and logn e = cases. In this figure, the y-ranges given by the 
measured ratios (thick lines) shift up and down as a result of decreases and increases in the 
abundances relative to oxygen, respectively; thus different parts of the theoretical curves 
would be in bold, moving the allowed temperature ranges (shown in the lower panel) to the 
left or right. Solar abundances here are taken to be the Cloudy 90 defaults (Ferland 1996; 
Grevasse & Anders 1989). 

In case (a), i.e. if the O Ylxa measurement of iVovi is correct, then the temperature 
and density of the absorber are tightly constrained in a metallicity-independent manner, 
and relative abundances for other elements can be estimated. Under this assumption, the 
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oxygen ion ratios are consistent within a range of n e = 10~ 4 - 7 — lCT 3 - 9 cm -3 ; however at 
solar abundances the Ne IX/O VII ratio demands higher temperatures than allowed by the 
O VI^/O VII ratio. Over this range of densities, the permitted abundances of neon, carbon, 
and nitrogen abundances relative to oxygen (that is, the range of abundances which produce 
line ratios consistent with both O Yl Ka /0 VII and O VIII/O VII) are then 0.6 <[Ne/0]< 2.2, 
-0.8 < [C/0]< 0.3, and [N/0]< 0.9. Note that supersolar [Ne/O] has also been observed 
in the z = absorber toward PKS2 155-304 (Nicastro et al. 2002). The improvement in 
the fit from supersolar [Ne/O] is shown in the right panel of Fig 6. Note, however, that 
in this case the discrepancy between the O Ylx a and 1032 A measurements becomes even 
more severe. Since the bulk of the O VI L v cannot be associated with the O VII due to the 
different Doppler parameters, the UV O VI component associated with the WHIM (hence 
the O VIko) must be substantially weaker than the Galactic O VI absorption; thus the 
discrepancy is correspondingly larger. 

On the other hand, if (b) the O VI#- a line is ignored then the relative abundances in 
the absorber can be adjusted such that the measured line ratios are consistent with either 
a low- or high-density absorber. As shown in Figure 6, a high-density, collisionally-ionized 
medium fits the data with Solar abundances. Assuming this is the case, the temperature 
is then completely constrained by the O VIII/O VII ratio at logT = 10 61-6 ' 2 K, and the 
relative abundances consistent with the O VIII/O VII ratio in this temperature range are 
-0.6 <[Ne/0]< 0.6, -0.6 <[C/0]< 0.3, and [N/0]< 0.4. 

The requirement that [Ne/O] be supersolar does not affect the viability of the partially- 
photoionized model: in both the ISM and low-z IGM, [Ne/O] is observed to be significantly 
larger than the solar value (e.g. Paerels et al. 2001; Nicastro et al. 2005a). This may be an 
intrinsic property of the enriched gas ejected into the IGM and ISM, or instead could be due 
to depletion of C, N, and O onto dust grains in supernova ejecta or quasar winds (Whittet 
1992; Elvis et al. 2002). If the dust destruction timescale is long and the IGM is continuously 
enriched by this latter mechanism, then the observed supersolar [Ne/O] would be expected. 
All of these enrichment scenarios are quite uncertain, but few (if any) are able to produce 
[Ne/0]< 0. Indeed, the solar neon abundance itself is quite uncertain since it is inferred from 
solar wind measurements. The increase in the solar neon abundance proposed by Drake & 
Testa (2005) is supported by these observations, and may provide another physical argument 
for the lack of subsolar [Ne/O]. 
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5.2. Where does the X-ray absorption originate? 



Assuming the VII absorption system is homogeneous, its radial extent can be esti- 
mated by calculating r ps iVf//(/m e ), where \i ~ 0.8, log(no/nH) = —3.13 (solar abundance), 
and 



The second term in the equation is approximately unity, since over the range of temperatures 
and densities implied from the O VIII/O VII ratio, O VII is the dominant ionization state 
by at least an order of magnitude; thus, logiV# = logiVovii + 3.13 — [O/H] = {20.37 — 
([O/H] + 1)} ± 0.11. The measurement error on iVovii is small compared to the uncertainty 
range in n e , so it can be disregarded. Assuming that (a) the O VI^ Q line does measure 
iVovi, the 2(T range of densities is —4.7 < logn e < —3.9, resulting in a radial extent of 
r = (0.8 — 4.9) x lcr([°/ H ] +1 ) Mpc. These radial extents are consistent with those expected 
from a Local Group medium or local filament interpretation for this absorption (Nicastro et 
al. 2002, 2003) and too large to be confined within a Galactic halo. The absorber extent could 
be marginally consistent with a Galactic corona if the metallicity is high (r = 80 — 490 kpc, 
2a limits at Solar metallicity); however, such a scenario seems implausible, particularly if 
this corona primarily consists of gas accreted from the metal-poor IGM. 

It is unlikely that this absorption system, if extragalactic, is spherically symmetric 
(particularly in a "local filament" interpretation). The total mass in the O VII system can 
be written as M tot — fx (4/3)7rr 3 (1.4nH?«H), where / (< 1) parameterizes the departure 
from spherical symmetry and lAnnmn ~ nnmu + nH e ^He- Replacing r with 7V#/(/m e ) and 
plugging in relevant values, 



This is several times larger than estimates of the total Local Group binding mass, e.g. 
M tot = (2.3 ± 0.6) x 10 12 M Q as calculated by Courteau & van den Bergh (1999). Baryonic 
matter should only contribute ~ 15% of this mass (if the baryon-to-dark matter ratio is 
equal to the cosmological value), so our estimate appears high. This discrepancy can be 
easily resolved with different values of [O/H] and /. For instance, if we assume an oxygen 
abundance of 0.3 times solar rather than 0.1, then the range of possible masses becomes 
2.0 x 1O U /M < M tot < 7.9 x 10 12 /M . Thus, unless [O/H] is very high or / < 1, the 
O VII absorber almost certainly accounts for a major fraction of the baryonic matter in the 
Local Group. This sightline may also be probing gas that is not gravitationally bound to the 
Local Group (i.e., the Local Filament), which may explain why our range of M tot extends 
to such high values. Although Collins et al. (2005) argue that a Local Group origin requires 






(3) 
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that the O VII absorbers contain too much mass, we see here that the total mass is in fact 
consistent with expectations. 

If this absorption does only trace the Local Group medium, then constraints on the 
extent of the absorber can be derived by assuming M tot = 0.15M LG ss 3.5 x 1O 11 M . In this 
case, r 3 f = 3M to t/(47rx 1.4/in e mH). Taking the 2a upper limit of logn H = —3.9, the 2a lower 
limit on Mlg = 11 x 1O 12 M , and assuming / = 1, we obtain a lower limit of r > 0.2 Mpc 
for the O VII absorption. Similarly, the upper limit on the radius is r < 0.6/ -1 / 3 Mpc. This 
seems somewhat small compared to the actual size of the Local Group, but once again is 
dependent on the geometry of the absorber. A value of / ~ 0.1 brings this upper limit 
more in line with the Mpc scales expected in the Local Group; this may indicate that the 
Mkn 421 line of sight probes an extended, filamentary WHIM distribution. This calculation 
also assumes that the density of the Local Group medium is constant with radius, when the 
actual density profile is more likely centrally concentrated. These measurements are also 
affected by the O VI discrepancy: in case (b), only the lower density limit of log(n e ) > —4.7 
(from the 2a O VI Ka upper limit) applies, so the upper radius and mass limits are still 
valid. Nevertheless, the consistency with the expected Local Group parameters is intriguing. 



5.3. Comparisons to Other Studies 

Kravtsov et al. (2002) used constrained simulations to study the properties of the Local 
Supercluster region; sky maps produced from this simulation (their Figure 5) show filamen- 
tary structures near the Mkn 421 direction, possibly corresponding to the observed absorp- 
tion. Additionally, they note that a Local Group medium would exhibit a low Doppler 
parameter (b < 60kms~ 1 ) out to distances of ~ 7 Mpc, consistent with our O VII measure- 
ment. 

The inferred properties of the X-ray absorption along this line of sight also appear similar 
to other observations: toward 3C 273, for example, Fang et al. (2003) find 5.36 < logT < 6.08 
and comparable O VII column density; however, their inferred O VII Doppler parameter is 
significantly higher than that toward Mkn 421: b > 100 km s _1 . On the other hand, the z = 
absorber toward PKS 2155 — 304 (Nicastro et al. 2002) exhibits a temperature consistent with 
the Mkn 421 absorber, yet inferred density about an order of magnitude lower. Compared to 
the two intervening filaments seen toward Mkn 421 (Nicastro et al. 2005a), the density of the 
z = absorption agrees with the derived lower limits (logn e > —5 for both filaments), but the 
filaments appear to exhibit higher temperatures than that derived for the local absorption. 
These variations along different lines of sight simply demonstrate the complex nature of 
the absorption, and the diversity of temperature and density environments produced in the 
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structure formation process (e.g., Kravtsov et al. 2002). 

The detection of X-ray absorption lines toward the Large Magellanic Cloud binary LMC 
X-3 by Wang et al. (2005) presents an important consideration for these results as well. The 
measured O VII and Ne IX equivalent widths and upper limits on O VIII and O VUr/s 
are all consistent with the same lines measured toward Mkn 421 (albeit with much larger 
statistical errors). Although this detection provides evidence of a hot intervening absorber 
between the Galaxy and LMC, it certainly does not rule out a primarily extragalactic origin 
for the Mkn 421 absorber. Any absorption, either Galactic or extragalactic, is likely to be 
inhomogeneous; thus, it's entirely plausible that the LMC X-3 sightline probes hot Galactic 
gas (perhaps enhanced by winds or outflows from both the Galaxy and LMC), while the 
absorption toward Mkn 421 is primarily due to low-density nearby WHIM gas. 



6. Summary and Future Work 

Through long-duration Chandra and FUSE observations of Mkn 421 in outburst, we 
have obtained unprecedented measurements of a variety of z ~ absorption lines, many of 
which likely arise in extragalactic, partially photoionized gas. A brief summary of our results 
follows. 

1. The relative strengths of the three O VII K-series lines imply 2a Doppler parameter 
constraints of 24 < b < 55kms _1 . This is inconsistent with the value of b = 80.6 ± 
4.2 km s" 1 derived for the Galactic low-velocity O VI, indicating that the O VIlv and 
O VII likely arise in different phases. The O VII b value is, however, consistent with 
the local IGM simulations of Kravtsov et al. (2002) out to distances of several Mpc. 

2. A weak high-velocity O VI 1032 A component also appears in the FUSE spectrum. 
Although its width is consistent with the O VII b measurement, the O VIhvc column 
density is too low to be associated with the O VII absorption unless T > 10 7 K (which 
itself is ruled out by the upper limit on b). The O VIhvc m &y thus represent a distinct 
third component along this line of sight. 

3. The column density inferred from the O VI K a line is a factor of ~ 4 higher than that 
measured from the O VI 1032 A transition. This may be due to unaccounted-for atomic 
physics effects, in which case the Ka line may provide a more accurate measurement 
of Aovi than the 1032 A line. We consider both cases: 

(a) If the O VI Ka line measures Aovi, then strong constraints on the tempera- 
ture, density, and relative abundances of the X-ray absorber can be derived: 
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T = 10 5 - 5 " 5 - 7 K and n e = 10" 47 - lO^cm" 3 , with allowed abundances of 
0.6 <[Ne/0]< 2.2, -0.8 <[C/0]< 0.3, and [N/0]< 0.9 (all 2a ranges). This 
range of densities, combined with iVovH) implies a total mass and extent consis- 
tent with those expected in the Local Group and/or Local Filament if the gas 
metallicity is low. However, in this case the O VI UV-X-ray discrepancy becomes 
worse since (due to the Doppler parameter constraints) only a small portion of 
the O VI L v line can be associated with the extragalactic X-ray lines. 

(b) If, instead, the O VI Ka line does not correctly measure iV vi, then the O VI 
associated with the O VII absorption must be fully blended with the Galactic 
1032 A O VI, and thus not measurable. In this case a lower density limit of 
logn e > —4.7 is found, which is consistent with either a Galactic or extragalactic 
medium. 

Much work remains to be done — both in order to better understand the data presented 
here, and to determine the true nature of the z ~ X-ray absorption. Higher signal to 
noise data along the Mkn 421 sightline would be useful to obtain better column density 
constraints, particularly on O VIII, O Vl^a, and Ne IX, and thus better constrain the effects 
of photoionization on the absorber. Data of comparable quality along other sightlines would 
be invaluable as well, both to probe other regions surrounding the Galaxy and to reconfirm 
the tantalizing O VI results presented herein. Higher-resolution simulations of the Local 
Group may allow us to determine the ionic column densities expected in the WHIM, and 
thus whether or not the observed absorption can possibly arise in the WHIM. Finally, more 
detailed modeling of O VI inner-shell transitions would shed a great deal of light on whether 
or not the O VI X-ray/UV discrepancy is real, and thus provide an invaluable framework for 
studying new (and existing) X-ray data. 
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Ne X Ka 
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VII K0 
VII K7 

Vlll^a 


24.737 
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13.447 
12.134 
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24.781 
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22.019 
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18.629 
17.768 
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19 1-1+0.03 
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1031.926 
1031.926 
1037.617 


1031.88 ±0.01 
1032.30 ±0.04 
1037.59 ±0.01 


152.3 ±7.0 

145.4 ± 11.6 


-13.7 ±2.6 
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-9.3 ±4.0 


270.9 ±7.9 
18.5 ±5.6 
164.9 ±8.4 


14.43 ±0.02 
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a From Verner et al. (1996), except O VI and O V which are from Schmidt et al. (2004) 

b Wavelength uncertainty is measured from fit or intrinsic LETG 0.02 A error, whichever is greater. For upper 
limits, the line position was frozen to the rest wavelength. 

c Error bars are la; upper limits are 2a. 

d Column densities are calculated using curve-of-growth analysis with b — 80.6 ± 4.2 km s -1 for C V, N VI, and 
O V, and the O VII Doppler parameter of b = 31 — 46kms~ 1 (la range) for all other lines. 

Note. — (1) The Ne X line lies within a detector feature, so only an upper limit on its equivalent width is 
given. (2) The O VI and O V A res t values are from laboratory measurements by Schmidt et al. (2004); theoretical 
oscillator strengths are taken from Pradhan et al. (2003). (3) There is an absorption line at — 760 kms -1 from the 
O V rest wavelength, but this is more likely O VII associated with Mkn 421. 




Fig. 1. — Portions of the Mkn 421 Chandra LETG spectrum (points) and the best-fitting 
model (histogram). Absorption lines at z ~ are labeled, and vertical tick marks indicate 
absorption from the z = 0.011 (solid) and z = 0.027 (dotted) intervening WHIM filaments 
(Nicastro et al. 2005a). 
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Fig. 2.— FUSE spectrum of Mkn 421 around the O VIAA1032, 1037 A absorption doublet 
(histogram). The dark curve shows the best-fit model with (solid line) and without (dotted 
line) the 1032 A HVC included. The inset plot shows the best-fitting Galactic (large Gaus- 
sian) and high-velocity (small Gaussian) components for the 1032 A line, plotted against 
velocity relative to the O VI rest wavelength. 
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Fig. 3. — Contours of allowed iVovn and b for the O VII K a (yellow), Kp (red), and 
lines. Shaded regions depict the la errors on W\, with the best-fit W\ line in the center 
of each region. The overlap between the K a and region provides a 2a lower limit of 
b > 24kms~ 1 , and the K a — Kp overlap implies b < 55kms _1 . The weighted average O VII 
column density of the three transitions over the la b range is log(iVovii) = 16.24 ± 0.11. 
The blue box depicts the la ranges in log(iVovii) and b. Also labeled on the top axis is 
logT max = log(m& 2 /2&;), the temperature inferred from a given Doppler parameter assuming 
purely thermal motion. 
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Fig. 4. — Contours of allowed iVovi and b for the O VI 1032 A (green), 1037 A (red), and 
putative Ka 22 A (yellow) lines. Contour (a) is derived from the 1032 A O VIlv equivalent 
width while (b) is from the measured FWHM; the intersection between the two green contours 
provides a tight constraint of b = 80.6 ± 4.2 kins -1 and log iVovi(cm~ 2 ) = 14.432 ± 0.016 for 
the O VIlv (shown as the blue cross). Note that the O VI Ka transition predicts Aovi about 
0.5 dex higher than the UV line, and this discrepancy cannot be explained by saturation. 
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Fig. 5. — Contours of constant abundance ratios for X-ray and UV oxygen absorption lines. 
Vertical bars denote the 2a range of temperatures inferred from the abundance ratio at 
each step in logn e . The horizontal dashed line is the 2a upper limit on the temperature of 
the O VII absorber from the Doppler parameter measurement. While the X-ray line ratios 
predict significantly different temperatures in the collisional ionization regime (logn e > —3), 
they are consistent with log7i e (cm~ 3 ) = —4.7 to —3.9 and logT(i^) ~ 5.5 — 5.7 (denoted 
by the black cross). Note that the O VIhvc/O VII ratio predicts a temperature far higher 
than the other derived limits, and thus the HVC is unlikely to be associated with the O VII 
absorption. 
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Fig. 6. — Same as Figure 5, but for ratios of several different ion abundances to VII: 
Yl Ka (green), O VIII (red), C VI (blue), Ne IX (cyan), and N VII (dotted black region). 
The black cross shows the range of logT and logn e derived from the O VIII/O VII and 
O VIa'q/O VII ratios. Assuming solar abundances (left panel), the O YIko an d Ne IX 
contours are inconsistent for logn e > —5, while all other ratios (except O VI^O VII) are 
consistent at logn e > —4.5. The right panel shows how a neon abundance shift of [Ne/0]= 1 
produces better agreement in the low-density regime. High-density models agree with the 
data only if the O Vl^a measurement is disregarded. 
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Fig. 7. — Models of ionic column density ratios for ions likely to arise in the Galactic ISM 
(assuming pure collisional ionization). Calculated ratios are shown as thin lines with ±2<r 
ratio measurements overplotted (thick segments). Upper limits are shown as dashed lines, 
and the temperature ranges derived from the different ion ratios are shown in the lower panel. 
The same-element ratios in this case provide strong constraints of logT > 5.4 (O V/O VI) 
and logT < 6.6 (N VII/N VI). 
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Fig. 8. — Same as Figure 7, for the ions likely originating in an extragalactic medium, with 
logn e = —3.9 (left) and logn e = (right). Solar abundances relative to oxygen are assumed, 
and all calculated ratios are relative to iVovn- The low-density case produces a consistent 
temperature solution of T ~ 10 5 7 K (from the overlap between the O VIII and O VIr- q 
bold regions), and [Ne/O] must be supersolar. The high-density case, on the other hand, 
provides a consistent solution with Solar abundances but the temperatures implied by the 
O VIII/O VII and O Vl Ka /0 VII ratios are inconsistent. 



